Pairwise entanglement in symmetric multi-qubit systems 
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The concurrence, a quantitative measure of the entanglement between a pair of particles, is de- 
termined for the case where the pair is extracted from a symmetric state of A*' two-level systems. 
Examples are given for both pure and mixed states of the A'^-particle system, and for a pair extracted 
from two ensembles with correlated collective spins. 



PACS numbers: 03.65.Ud, 03.67.Lx, 75.10.Jm. 
I. INTRODUCTION 

Various proposals exist for the preparation of multi- 
particle entangled states, and a number of these states 
have been pointed out to be particularly easy to prepare 
and to have special and useful properties. Since entan- 
glement is defined as a property of the whole ensemble 
of particles, it is not immediately clear whether a sub- 
ensemble of particles, drawn at random from the original 
ensemble will also be in an entangled state, or whether 
the trace over some particles will destroy the quantum 
correlations in the system. In this paper we consider 
the simple questions whether a random pair of particles, 
extracted from a symmetric state of N two-level systems 
will be in an entangled state or not, where by symmetric, 
we assume symmetry under any permutation of the par- 
ticles. Entangled states constitute a valuable resource 
in quantum information processing j^], and the trans- 
fer of entanglement between few qubits and the quasi- 
continuous variables by which we describe many-particle 
systems, may become an important ingredient in, e.g., 
quantum data-storage and inter-species teleportation. 

The paper is organized as follows. In Sec. II, we 
present the concurrence, introduced by Wootters [H,D, 
who demonstrated its one-to-one correspondence with 
the entanglement of formation of a pair of qubits. In 
Sec. Ill, we show how the density matrix of a pair of 
qubits can be expressed in terms of expectation values 
of collective spin operators on the multi-qubit state. In 
Sec. IV, we analyze three examples of pure states of 
the N particles: spin coherent states, Dicke states, and 
spin squeezed states. In Sec. V, we consider an example 
of a mixed state with thermal entanglement Q-^, and 
we show examples where the pairwise entanglement de- 
pends on the temperature of the system. Finally in Sec. 
VI, we assume two separate ensembles in an Einstein- 
Podolsky-Rosen state of correlated angular momentum 
components, and we show that a single pair with an atom 
from each ensemble will be in an entangled state. 



II. TWO-PARTICLE DENSITY MATRICES AND 
ENTANGLEMENT 

It is easy to check if a pure state of two quantum sys- 
tems is an entangled state or not by simply observing 



the eigenvalues of the reduced density matrix of either 
system. It is also possible to quantity the amount or de- 
gree of entanglement of the state ||], _B = — log2 r^, 
which presents the asymptotic ratio between n and m, 
where n is the number of pairs in the desired state, syn- 
thesized from m pairs of maximally entangled states. 

For a mixed state with density matrix /O12, a simi- 
lar measure can be defined as the minium value of the 
weighted average of E over wave functions by which the 
two-particle density matrix can be written as a weighted 
sum. It is necessary to search for the minimum, since 
P12 can be written in many ways as a weighted sum 
of pure state projections. In the general case, this is a 
highly non-trivial task, as is the determination whether 
the state is entangled at all. For two qubits, however, 
entanglement is equivalent with the non-positivity of the 
partially transposed density matrix ]lO| , and the entropy 
of formation can, magically, be obtained as a simple an- 
alytical expression 



E = h 



l + Vl^ 



(1) 



where h{x) — — a;log2 x ^ {1 — x) log2(l — x), and where 
the concurrence, C, is defined as 



C = max {0, Ai — A2 — A3 — A4} 



(2) 



where the quantities A^ are the square roots of the eigen- 
values in descending order of the matrix product 



£'12 = Pl2{<yiy ® ^2y)p*i2{^'^V ® ^^y)- 



(3) 



In 



^ p'{2 denotes the complex conjugate of pi2, and 
ly are Pauli matrices for the two-level systems. The 



eigenvalues of Q12 are real and non-negative even though 
Q12 is not necessarily Hermitian, and the values of the 
concurrence range from zero for an unentangled state to 
unity for a maximally entangled state. 



III. DENSITY MATRIX FOR A PAIR OF QUBITS 
FROM A MULTI-QUBIT STATE 

A two-qubit reduced density matrix which is symmet- 
ric under exchange of the two systems can be written 
as 
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(4) 



where the matrix 

elements in the basis {|00), |01), |10), |11)} can be rep- 
resented by expectation values of Pauli spin matrices of 
the two systems 

w± = ^ (1 ± 2(criz) + (cri^a-2z)) , 

X± = + ) ± (cti + CT2z)), 



- (1 - (ai^a2z)) 



w 

y = ^o-i+o'2 



{o'lyO-2y) + i2{(Jix<72y)) ■ 



(5) 



We now consider the entanglement of two qubits ex- 
tracted from a symmetric multi-qubit states. If only sym- 
metric qubit states are considered, we can describe the 
state of the A^-qubit system in terms of the orthonormal 
basis \S,M){M = -S,-S +1, ...,S) with 5 = 7V/2. The 
states \S^M) are the usual symmetric Dicke state [0, 
i.e., eigenstates of the collective spin operators S'^ and 
S^, defined as 



1 ^ 

2 ^ 



x,y,z. 



(6) 



N{N - 1) ' 

{[S+,Sz]+) 
N{N- 1) ' 



(10) 



where = AB + BA is the anticommutator for 

operators A and B. 

From Eqs.(|) and (|lO|) we may thus express the density 
matrix elements of pi2 in terms of the expectation values 
of the collective operators. 



v± 
x± 
w 



-2N + 4(52) ± A{Sz)iN - 1) 
4A^(A^- 1) ' 
iN-l){S+)±{[S+,SzU) 
2N{N-1) 
- 4(5f ) 
4:N{N-1) ' 
2(5g + S^y) - N 
2N{N~1) ' 
(gg-5g)+»([5.,5,]+) 
N{N - 1) 



N{N- 1) 



(11) 



IV. PURE MULTIQUBIT STATES 

In this section we study three examples, where the N 
two-level systems are described by a pure state which is 
invariant under permutation of the particles. 



For later use it is convenient to define the number op- 
erator Af — Sz + N/2 and number states as 



\n)N = \N/2,~N/2 + n)N, 
Af\n)N = n\n)N- 



(7) 



The eigenvalue n of the number operator A/" is the number 
of qubits in the state |0). For example, the states |0)Ar 
and |l)Ar are explicitely written as 



|0)^ 

|1>A^ 



I111...1), 

^ =(|011...1> 



I111...0)). 



N 



1101. .1) 



(8) 



(9) 



|l)Ar is also called an iV-qubit W state |12|,p|. 

Due to the symmetry of the state under exchange of 
particles we have 



(o-ic 



2{S^) 

N 
{S+) 

N ' 
4(^^) 



N 



N{N - 1) 



A. Spin coherent states 

The spin coherent state |l^ is obtained by a rotation of 
the spin state \S,M = S), which in turn is the product 
state of all TV particles in the |0) state. Hence it is a 
separable state. It is still interesting to go through the 
above procedure and to insert the explicit expression of 
the spin coherent state [|l3|. 



h) = (i + H') 



N 

2\~N/2 

n=0 



1/2 



??"|f^)iV, 



(12) 



where rj is chosen real in the following. By a straightfor- 
ward calculation from Eqs.(|Tl|) and (p^, we find 



1 



Pl2 



(l+??2)2 



/ f]'^ rf rj^ 

rj rj 1] rj 

rj rj rj rj 

rj^ rj rj 



(13) 



1 



which is in agreement with our observation that the two- 
particle state is really a product state of two rotated spin- 
i particles in the states {rj\Q) / y/l + rf . The matrix 
product Qi2 is found to be a 4x4 matrix of zero's, reveal- 
ing the role of the ay Pauh matrices in (||): pi2 is the 



2 



projection operator on spin states with a definite direc- 
tion in the xz-plane, the application of ay is equivalent 
to a 180° rotation in the xz-plane, and gi2 is therefore 
the vanishing product of projection operators on two or- 
thogonal subspaces. Naturally, the concurrence vanishes 
in this case, C — 0; there is no pairwise entanglement in 
the spin coherent state. 



B. Dicke State \N/2, M) 

The Dicke states, defined as effective number states 
above, are states with a definite number of particles oc- 
cupying the internal states |0) and 1 1) . Such states may in 
principle be prepared in an atomic physics experiment by 
Quantum Non-Demolition detection of the atomic popu- 
lations by phase contrast imaging of the atomic sample 
1^,^. By rotation of all spins, a separable spin coherent 
state is first prepared with a binomial distribution on the 
various Dicke states, cf., Eq.jl^), and experiments have 
already demonstrated a factor 3 reduction in the variance 
of the populations after such a detection |Q . 

From Eq.(|ll|), it is easy to see that the reduced density 
matrix pi2 is given by 
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with matrix elements 

{N ±2M)(N -2±2M) 



(14) 



v± 



AN{N - 1) 



AN{N - I)' 



(15) 



The concurrence of a simple density matrix of the form 
( ^ ) with x±, u = 0, and y ^ y* is given by [T^ 



C ~ 2 max{0, y — 



(16) 



where we have used that 2y + w+ + w- = 1 to pi2 . Now 
substituting Eq.(|l5|) into ([l^), wc explicitly obtain 



C 



1 



-{N^ ~ AM- 



/r2 



2N{N- 1) 

-^/{N^ - AAP)[{N - 2)2 - AAP]}. (17) 



The values of C for different N and M are illustrated in 
Fig. 1. For any Dicke state except the ones with maxi- 
mum \M\, if one extracts two particles, they will be in an 
entangled state. We also observe that the concurrence is 
nearly a constant in the neighborhood of M = 0. 
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FIG. 1. The concurrence in the Dicke state for diflterent 
number N. N = 15(open circles), N — 25 (crosses), and 
TV = 30 (open square). 

The variation of Af around M = is small for an 
initial binomial distribution with this mean value, and 
the concurrence will be very close to the exact result, 
C — 1/{N — 1) for M = 0, irrespective of the outcome of 
a QND measurement of M. 

The Dicke states |7V/2, M = ±{N/2 - 1)) have a con- 
currence of C = 2/N. These states are identical with the 
W state (see Eq.(||)), which are known to be the sym- 
metric states with the highest possible concurrence [18|. 



C. Kitagawa-Ueda state 

In 1993, Kitagawa and Ueda proposed a nonlinear 
Hamiltonian x^x i^i order to generate spin squeezed 
states [ p^ . This effective Hamiltonian may be realized in 
ion traps [ pO| , where it was already implemented in order 
to produce multi-particle entangled states (of four parti- 
cles) 1 21 , and it may be implemented in two-component 
Bose-Einstein condensates as a direct consequence of the 
coUisional interactions between the particles |Q , see also 

il- 

When the Hamiltonian H — x^l is applied to the 
many-particle system, which has been prepared in the 
product state |0)Ar = |111,...,1), the wave function at 
time t is obtained as 



l*(t)> 



|0) 



N- 



(18) 



Using the results obtained in ||l9[ the following expecta- 
tion values are obtained (/i = 2xt) 



(Sx) = (Sy) = 

(Si) = N/A 



N- 



HI) 



3 



(S^) = - {N^ + N~ N{N - 1) ( 



{SD = -{N^ + N + N{N - 1) 

([5+,^,]+) =0 

{[S., Sy]+) = In{N - 1) ^ sin ^. (19) 

We are now able to determine the two-particle density 
matrix, which is on the form 



Pl2 



/ v+ u* 

w w 

w w 

\u w_ 



(20) 



with matrix elements given by Eq. ( |ll| ) . The combination 
of Eqs.(pl]) and ( p^ ) gives explicitly the matrix elements. 

From Eqs-d) and (^), the concurrence for the matrix 
( pO| ) is obtained as 



C = 



2 max(0, |u| — w), if 2w < y/v^v^ + \u\; 
2 max(0, w — y/v-^-V-), if 2w > y/v+v^ + \u\. 



(21) 



The concurrence of the spin squeezed states is given by 
analytical expressions in the argument fj, = 2xt, which 
are too lengthy to present her. In Fig. 2 we present the 
results numerically: If two atoms are extracted at ran- 
dom from spin squeezed samples they will be in a mutu- 
ally entangled state. We observe that the concurrence is 
symmetric with respect to n = tt. At this special point 
of = 77 the A^-particle GHZ state is produced |20|, and 
it has no pairwise entanglement. 




FIG. 2. The concurrence as a function of n for different 
number N. N = 3(solid line), iV = 4 (dashed line), iV = 5 
(open circles), A'^ = 6 (open diamonds), and N = 7 (crosses). 



V. MIXED MULTIQUBIT STATES AND 
THERMAL ENTANGLEMENT 

An interesting and novel type of thermal entangle- 
ment was introduced and analyzed within the Heisen- 
berg XXX 0, XX g, and XX Z @ models as well 
as within the Ising model in a magnetic field The 
state of the system at thermal equilibrium is represented 
by the density operator p(T) = exp (—H/kT) /Z, where 
Z =tr[exp {—H/kT)] is the partition function, H the sys- 
tem Hamiltonian, k is Boltzmann's constant which we 
henceforth take equal to unity, and T the temperature. 
As p(T) represents a thermal state, the entanglement in 
the state is called thermal entanglement Unlike in 
standard statistical physics where all properties are ob- 
tained from the partition function, determined by the 
eigenvalues of the system, entanglement properties re- 
quire in addition knowledge of the eigenstates. The ana- 
lytical results in the previous studies on thermal entan- 
glement are only available for two and three qubits 
g. Here we consider pairwise entanglement in the mul- 
tiqubit systems. 



A. Isotropic Heisenberg model 

We consider the A^-qubit isotropic Heisenberg Hamil- 
tonian 



N 



(22) 



The positive (negative) J corresponds to the antiferro- 
magnetic (ferromagnetic) case. In this model all particles 
interact with each other. 

By using the collective spin operators, the Hamiltonian 
is rewritten as 



H = j{Sl + Sl + Si) = JS^ 



(23) 



up to a trivial constant. 

Unlike pure states, the symmetric multi-particle den- 
sity matrix does not only populate the fully symmetric 
Dicke states, and we have to determine the number of 
collective spin-S" states for each S. Write S as (A/2 — k), 
we know that for A: = 0, a single irreducible representa- 
tion exists: the A -I- 1 fully symmetric Dicke states with 
S — N/2. The number of irreducible representation with 
S — N /2 — 1 is obtained by noting that their maximum 

N 



M value is also N/2~ 1, and a total of 



1 



N states 



exist with precisely one particle in the |1) state. One of 
these belong to the S — N/2 irreducible representation, 
and the remaining A^ — 1 states must have S = N /2 — 1. 
This argument can now be repeated to obtain the num- 
ber of states with S = A/2-2 and M = N/2-2, i.e., the 
number of 5' = A^/2 — 2 irreducible representation, etc., 
until all 2^ states of the system have been accounted for. 
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The isotropic Hamiltonian only depends on 5^, and 
knowing the multipHcity of each value of this quantity 
we write the partition function 



N/2 

Z = ViVfc[2(7V/2-fc) 



fe=0 



-l3,J(N/2-k)(N/2-k+l) 



(24) 



where Nk — 



N 
k-1 
f N 

argument. We assume I _^ 

The reduced density matrix for two qubits is 



follows from the above 



0. 



Pl2 



/ w 

Q w y Q 

y w 

V i; 



(25) 



The matrix ele- 



with matrix elements given by Eq.(|l] 
ment v = v± since (Sz) = 0. 

Due to the symmetry property {S^) — (Sy) — {Si), we 
only need to determine 

N/2 N-2k 

{SD^Y.^^ E + (26) 

k=0 m=0 
^-pj(N/2-k)(N/2^k+l) 

And from Eq. ([T^) , the concurrence is obtained as 
1 



— max{0, 
■ max{0, 



2A^(A^~ 1) 

2\l{Sl + 5^) -N\-N'^ - 4(52)} 
1 



2iV(iV- 1) 

2|4(5f) - iV| - TV^ _^ 2iV - 4(51)}. 



(27) 



To identify the sign oi A = 2|4(S'2) - iV| - TV^ + 2A^ - 
4(5^) in (H^), we consider the case where 4(5"^) > for 



which A = 4(52) - 7V2. Since (5^) < if (f + l) < i^, 
we always have 4(5^) — N"^ < 0, and there is no pairwise 
entanglement. In the opposite case where 4(5^) < iV, we 
have A = 4iV - 12(5^) - TV^ since (5^) > 0, we have 
A < if iV > 4. For case of iV = 3, we have shown that 
the pairwise thermal entanglement is absent from both 
the antiferromagnetic and ferromagnetic isotropic model 

So we conclude that there is no thermal entanglement 
for TV > 3 in the isotropic Heisenberg model. The case of 

= 2 is discussed in detail in Ref. |j] and it is shown that 
there is no thermal entanglement for the ferromagnetic 
case. In order to observe the pairwise entanglement in 
the multiqubit system, now we consider the anisotropic 
Heisenberg model. 



B. Anisotropic Heisenberg model 

The anisotropic Heisenberg Hamiltonian is given by 
Ha = J [si + 52 + A52) = J52 + J(A - 1)52, (28) 

where A is the anisotropy parameter. Obviously the 
Hamiltonian Ha reduces to Hi when A = 1, and Ha 
yields the XX model when A = 0. 

The concurrence is still given by (p7|), but the partition 
function and the relevant expectation values now become 

N/2 N-2k 

k=0 m=Q 

-I3J(N /2-k)(N /2-k+l) 



-p,J{A-l)(m-N/2+kY 



(29) 



{Si 



xe 

N/2 



k=0 



N-2k 

E 

m=0 



{Si 



Si) = 



-f3J{A-l}(m-N/2+kf 
-pJ(N/2-k){N/2-k+l) 1^ 

[(A/2-fc)(Ar/2-fc + l) 



N-2k 



xe 

N/2 
k=0 

-{m- N/2 + kf 



(30) 



m=0 



xe 



-0J{A-l){m-N/2+ky 
-0J{N/2-k){N/2-k+l) 



/z. 



(31) 



This model leads to pairwise entanglement, as shown 
by the numerical results presented in Figure 3 as func- 
tions of the reciprocal temperature, x = fiJ . For N = 2 
we observe that the concurrence is symmetric with re- 
spect to X = 0, which is consistent with the result in Ref. 
g. In other words, the thermal entanglement appears 
for both the antiferromagnetic and ferromagnetic cases. 
However for A^ > 3, the thermal entanglement only exists 
for the ferromagnetic case. We observe a critical value 
of X, after which the entanglement vanishes. And the 
critical value increases as A^ increases. 
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FIG. 3. The concurrence as a function of a; = /3J for dif- 
ferent number A'' in the XX model (A = 0): A'^ = 2(open 
circle), N — 5 (crosses), A*' = 15 (solid line), and N = 25 
(dashed line). 

Within the above framework we may also consider 
more general models such as 



Hg = J S'^ 



(32) 



where f{Sz) is an arbitrary analytical function of Sz- As 
the operator f{Sz) commutes with 5*^, similar analytical 
results for the concurrence can be obtained and the ther- 
mal entanglement can be readily generated for special 
choices of f{Sz)- 



VI. EPR-CORRELATED ENSEMBLES 

Finally we consider two EPR-correlated ensembles. 
This state is not invariant under any permutation of par- 
ticles, but only under those permutations that exchange 
particles within each ensemble, and it is furthermore 
characterized by the correlations between the samples 1 
and 2: 



(Ji. - J2x)|*) = 0, 
(Jl,+J2,)|*> =0. 



(33) 
(34) 



A state that obeys Eqs.(p3|) and ( |34| ) can in principle be 
obtained by successive QND detection of the observables 



Jix ~ J2x and Jiy + J; 



J2y 



25|. Equivalently the above 



equations can be written as 

{Jl+ - J2-) 
{Jl- - J2+) 



Vf) =0, 

^) = 0. 



(35) 
(36) 



It is easy to check that a solution of the above equation 
is the EPR-correlated state 



1 ^ 

\n)N ® \n)N 



(37) 



V= -{l±2{(Tiz) + {aiz(J2z)), 

^ n / w ^ {•hzJ2z) 



4^- 4 ^2 ' 

^ (J1+J2+) 

The concurrence is given by 

C ^ 2max{0, \u\ - w} 

(J1+J2+) -I- {Jiz'hz) 



(39) 



2 max{0, ■ 



(40) 



The expectation values of Ji+ J2+ and JizJ2z are readily 
obtained in the state j^"), and we find that C — 1/N. The 
pair of particles is in an entangled state. If the ensembles 
are really macroscopic, as in the entanglement is, 
however, very weak. 



VII. CONCLUSIONS 

The purpose of this paper has been to point out that 
multi-particle entanglement quite typically implies pair- 
wise entanglement within the sample. We showed that 
the two-particle density matrix is readily expressed in 
terms of expectation values of collective operators, in the 
case of symmetrical states of the many-particle system, 
and we provided the value of the concurrence for a num- 
ber of examples. These results confirmed and general- 
ized results obtained, e.g., on the pairwise entanglement 
in systems with definite (TV = 3,4) numbers of particles. 

The entropy of formation, and the very issue of entan- 
glement, is highly non-trivial for situations dealing with 
more than two particles, and for mixed states of systems 
with dimensions higher than 2. Studying and optimiz- 
ing the two-particle concurrence in systems with many 
particles may be a useful way to learn about the more 
complicated case. 



And it also satisfies {Jiz ~ J2z)\'^) — 0. 

Now we consider the entanglement of two qubits, which 
belong to different ensembles, we first identify the two- 
qubit reduced density matrix : 



P12 



f V Q u* 

w 

w 

u Q V 



(38) 
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